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Abstract 

Dynamical linked cluster expansions are linked cluster expansions with 
hopping parameter terms endowed with their own dynamics. This amounts to 
a generalization from 2-point to point-link-point interactions. An associated 
graph theory with a generalized notion of connectivity is reviewed. We discuss 
physical applications to disordered systems, in particular to spin glasses, such 
as the bond-diluted Ising model and the Sherrington-Kirkpatrick spin glass. 
We derive the rules and identify the full set of graphs that contribute to the 
series in the quenched case. This way it becomes possible to avoid the vague 
extrapolation from positive integer n to n = 0, that usually goes along with 
an application of the replica trick. 
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1 Introduction 



Linked cluster expansions (LCEs) have a long tradition in statistical physics. Orig- 
inally applied to classical fluids, later to magnetic systems ([0,@,@ and references 
therein), they were generalized to applications in particle physics in the eighties [|J. 
There they have been used to study the continuum limit of a lattice $ 4 field theory 
in 4 dimensions at zero temperature. In [| they were further generalized to field 
theories at finite temperature, simultaneously the highest order in the expansion 
parameter was increased to 18. Usually the analytic expansions are obtained as 
graphical expansions. Because of the progress in computer facilities and the devel- 
opment of efficient algorithms for generating the graphs, it is nowadays possible to 
handle of the order of billions of graphs. The whole range from high temperatures 
down to the critical region becomes available, and thermodynamic quantities like 
critical indices and critical temperatures are determined with high precision (the 
precision is comparable or even better than in corresponding high quality Monte 
Carlo results) 0-0. An extension of LCEs to a finite volume in combination with 
a high order in the expansion parameter turned out to be a particularly powerful 
tool for investigating the phase structure of systems with first and second order 
transitions by means of a finite size scaling analysis [[10 . 



Linked cluster expansions are series expansions of the free energy and connected 
correlation functions about an ultralocal, decoupled theory in terms of a hopping 
parameter K. The corresponding graphical representation is a sum in terms of 
connected graphs. The value of K parametrizes the strength of interactions between 
fields at different lattice sites. Usually they are chosen as nearest neighbours. In 
contrast to the ultralocal terms of a generic interaction we will sometimes refer to 
hopping terms as non-ultralocal. 

In this paper we develop dynamical linked cluster expansions (DLCEs). These 
are linked cluster expansions with hopping parameter terms that are endowed with 
their own dynamics. Such systems are realized in spin glasses with (fast) spins and 
(slow) interactions [|TTJ-||13[. They also occur in variational estimates for SU(N)- 
gauge-Higgs systems, cf. |]14j| . Like LCEs they are expected to converge for a large 
class of interactions. 

Formally DLCEs amount to a generalization of an expansion scheme from 2-point 
to point-link-point-interactions. These are interactions between fields associated 
with two points and with one pair of points called link. The points and links 
are not necessarily embedded on a lattice, and the links need not be restricted 
to nearest neighbours. We have developed a new multiple-line graph theory in 
which a generalized notion of connectivity plays a central role. Standard notions of 
equivalence classes of graphs like 1-line irreducible and 1-vertex irreducible graphs 
have been generalized, and new notions like 1-multiple-line irreducible graphs were 
defined in order to give a systematic classification. 

The paper is organized as follows. In Sec. 2 we specify the models that admit a 
DLCE. We introduce multiple-line graphs and explain the idea behind the abstract 
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notions of multiple-line graph theory. Detailed definitions of multiple-line graphs , 
related notions and the computation of weights are given in Sect. 3. Sect. 4 treats 
the issue of renormalization in the sense of suitable resummations of graphs. Ap- 
plications to spin glasses are presented in Sect. 5. There it is of particular interest 
that DLCEs allow for the possibility of avoiding the replica trick. In the quenched 
limit we derive that DLCEs must be restricted to a subclass of the correspond- 
ing graphical expansion, so-called quenched DLCEs (QDLCEs). We also list some 
examples for models whose phase structure is accessible to QDLCEs. To these be- 
longs in particular the bond diluted Ising model. Sect. 6 contains the summary and 
conclusions. 



2 A Short Primer to DLCEs 

In this section we first specify the class of models for which we develop dynamical 
linked cluster expansions. Next we illustrate some basic notions of multiple-line 
graph theory, in particular the need for a new notion of connectivity. 

By A we denote a finite or infinite set of points. One of its realizations is a 
hypercubic lattice in D dimensions, infinite or finite in some directions with the 
topology of a torus. Ai denotes the set of unordered pairs (x, y) of sites x, y G A , 
x 7^ y, also called unoriented links, and A x a subset of A^ 

We consider physical systems with a partition function of the generic form 

Z(H,I,v) = ex.pW(H,I,v) 

= M j P0PC/exp(-5(0,C/,^))e X p(J]//(x)0(x) + J]/(OC/(O), (1) 

zeAo ieM 

with measures 

V(f) = Y[ d<t>{x) , VU = Y[ dU(l) (2) 

and action 

S(<l>,U,v) = J2 S °(<f>(x)) + J2 Sl ( U ( l »-l E v{x,y)<t>{x)U{x,y)4>{y), (3) 

with non-ultralocal couplings 

v(x,y) = v(y,x) ^ only for (x, y) G Ax, 

in particular v (x, x) = 0. (4) 

For later convenience the normalization via M is chosen such that W[0, 0, 0] = 0. 

The field <p(x) is associated with the sites x G A and the field U (/) lives on the 
links I E Ax, and we write U(x,y) = U(l) for / = (x,y). For definiteness and for 
simplicity of the notation here we assume (f)(x) G R and U{1) G R. In our actual 
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applications to spin glasses the 0s are the (fast) Ising spins and the Us G R are the 
(slow) interactions. The action is split into two ultralocal parts, S° depending on 
fields on single sites, and S l depending on fields on single links I G Ai. For simplicity 
we choose S 1 as the same function for all links I G Ai . We may identify Ai with the 
support of v, 

Ai = {l = {x,y)\v{x,y)^0}. (5) 

The support of v(x,y) need not be restricted to nearest neighbours, also the precise 
form of S° and S 1 does not matter for the generic description of DLCEs, S° and 
S 1 can be any polynomials in and U, respectively. The only restriction is the 
existence of the partition function. 

Note that the interaction term v(x,y) <f)(x) U(x,y) <p(y) contains a point-link- 
point-interaction and generalizes the 2-point-interactions v(x,y) <f)(x) 4>{y) of usual 
hopping parameter expansions. The effective coupling of the fields has its own 
dynamics governed by S X (U), the reason why we have called our new expansion 
scheme dynamical LCE. 

Dynamical linked cluster expansions are induced from a Taylor expansion of 
W(H, I, v) = In Z(H, I, v) about v = 0, the limit of a completely decoupled system. 
We want to express the series for W in terms of connected graphs. Let us consider 
the generating equation 

dW/dv(xy) = 1/2 < <f>(x)U(x,y)(f){y) > 

= 1/2 (W H y)H(y) + W H(x)Hiy) W I( 

x ,y) 

+ WH(x)I(x,y)WH(y) + Wi( XtV )H(y)WH(x) 

+ W H ( x )W H (y)Wi (x>y) y (6) 

Here < • > denotes the normalized expectation value w.r.t. the partition function 
of Eq. (|I|). Subscripts H(x) and I(x,y) = I(y,x) = 1(1) denote the derivatives of 
W w.r.t. H(x) and I(x,y), respectively. 

Next we would like to represent the right hand side of Eq. (Q) in terms of con- 
nected graphs. Once we have such a representation for the first derivative of W 
w.r.t. v, grapical expansions for the higher derivatives can be traced back to the 
first one. 

For each W in Eq. @ we draw a shaded bubble, for each derivative w.r.t. H 
a solid line, called a 0-line, with endpoint vertex x, and for each derivative w.r.t. 
I a dashed line, called a ?7-line, with link label I = (x,y). The main graphical 
constituents are shown in Fig. 1. Two 0-lines with endpoints x and y are then 
joined by means of a dashed ?7-line with label /, if the link I has x and y as its 
endpoints, i.e. I = (x,y). According to these rules Eq. (^), multiplied by v(x,y) and 
summed over x and y, is represented by Fig. 2. Note that, because of the Taylor 
operation, each solid line from x to y carries a factor v(x,y). 

Since the actual need for a new type of connectivity is not quite obvious from 
Fig. 2, because Eq. ffl) does not contain higher than first order derivatives w.r.t. /, 
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let us consider a term 




Figure 1: Graphical representation of the derivatives of W(H,I,v). (a) n-point 
function d n W/dH( Xl ) ■ ■ ■ dH(x n ), (b) n-link function d n W/dI(h) ■ ■ ■ dI(Q. 




Figure 2: Generating equation of the graphical expansion of DLCEs. The solid line 
in each graph carries a propagator v(x, y). A dashed [/-line with label / intersects a 
solid line with endpoints x and y if I — (x,y). 

W H{x) W H{y) W H{r) W H{s) Wj 

(x,y)I(r,s) (7) 

occurring in the second derivative of W w.r.t. v(x,y), v(r,s). According to the 
above rules this term would be represented as shown in Fig. 3a. While the 2 
vertices in the last term of Fig. 2 are connected in the usual sense via a common 
(solid) line (the dashed line with an attached bubble could be omitted in this case), 
the graph in Fig. 3a would be disconnected in the old sense, since neither x nor y 
are line-connected with r and s, but -as a new feature of DLCE graphs- the lines 
from x to y and from r to s are connected via the dashed lines emerging from a 
common bubble shown in the middle of the graph. As we see from Fig. 3a, we 
need an additional notion of connectivity referring to the possibility of multiple-line 
connectivity. While the analytic expression is fixed, it is a matter of convenience 
to further simplify the graphical notation of Fig. 3a at v — 0. Two possibilities are 
shown in Fig. 3b and Fig. 3c. To Fig. 3b we later refer in the formal definition 
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of the new type of multiple-line connectivity. In the familiar standard notion of 
connectivity two vertices of a graph are connected via lines. The vertices are line- 
connected. Already there, in a dual language, one could call two lines connected 
via vertices. The second formulation is just appropriate for our need to define when 
two lines are connected. The corresponding vertices mediating the connectivity of 
lines are visualized by tubes, in Fig. 3b we have just one of them. The tubes should 
be distinguished from the former type of vertices represented as full dots which are 
connected via bare 0-lines. In Fig. 3c we show a simplified representation of Fig. 3b 
that we actually use in graphical expansions. 




Figure 3: Representation of WH(x)WH(y)WH( r )WH( s )Wi( X} y)i( r}S y (a) according to 
the rules of Fig. 1 and 2, (b) same as (a), but at v — and simplified for a formal 
definition of multiple-line connectivity, cf. section 3, (c) same as (b), but for use in 
the actual graphical representations. 

The derivative terms have to be evaluated at v = 0. For v = we have a 
decomposition of W according to 

W(H, I,v = 0)=J2 W°{H{x)) + (8) 



with 



and 



f°° dcfiexp (-S°(4>) + H(j)) 
J_ oo d0exp(-S' (</>)) 



ex P W\l) = Z\I) = J -°° Pl \ (10) 

In Eq.s (^],|TUp we have omitted any single site or single link dependence, because 
we assume that S° and S l are the same for all x G Aq and i G Aj, respectively. 
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Therefore, at v = 0, the only non- vanishing derivatives of W are 
W I d n W°(H( Xl )) 

W H (x 1 )H(x2)...H(x n )\ v=0 ~ QH(xi) n " { > 

and 

vVl{h)i{h)...l{i m )\ v=0 ~ — dl(h) m dh ' 1 ^- 1 ^ { > 

but mixed derivatives w.r.t. H and I vanish. As anticipated in Fig.s 3b and 3c, for 
v = we replace the dashed bubbles and graphically distinguish between bubbles 
with Klines and bubbles with [/-lines. We define 

oc ( d"W°(H) \ 

- = «• = lrMF-J M (13) 

for a connected n-point vertex with n > 1 bare 0-lines emerging from it and 





d v W\I) 



" = = ( ^rf 1 ) ( 14 ) 

/=0 



for a connected z/-line consisting of u bare lines. If i/ = 1, we often omit the dashed 
line. If the bare lines of a z/-line are internal 0-lines, they get vertices attached to 
their endpoints, if they are external [/-lines, no vertices will be attached. 

Let V denote the lattice volume in D dimensions. The Taylor expansion of W 
about v = to second order in v then reads 

W(H,I,v) = W(H,I,v = 0) 



+ 



v{x,y)-W H{x) W H{y) W I 



x,y£A 



+ \ Yl \v{x,y)v(r,s) 



4 Wn(y) Wh( s ) W H (r)H(x) Wl{x,y) Wlfag) 
+ 2 WH(x)H(r)WH(y)H(s)Wi( Xi y)Wi( r>s ) 
+ 4 W H ( y ) W H ( S ) W H (r)H(x) Wi( Xj y)i(r l8 ) 
+ 2 WH(r)H(x)W H (y)H{s)Wi( rt y)i( XiS ) 

+ Wh(x) Wh{ v ) Wh(t) Wh( s ) Wi( Xt y)i( rjS ) ) 

+ 0(v 3 ), (15) 
where we have used that v (x, x) = 0. For each W in the products of Ws we now 



insert Eq.s (0) , (|12|) . 
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If we choose v in a standard way as next-neighbour couplings 



D-l 



J x-fi,y) 



(16) 



v(x, y) = 2K22(6 x+ {i t y + 6 X 

with fi denoting the unit vector in /i-direction, Eq. fljjD becomes in a graphical 
representation at H = I = 



W{0,0,v) 
V 



= (2K) l - {2D) * 
+ {2Kf [\ (2Df 



+ 2 (2^) 2 
+ 1 2(2D) 




+ 4 ( 2Z) ) 



+ 4 (2D) 




(17) 



+ 0(K 3 ). 

For clarity, here we have written explicitly the topological symmetry factors and 
the lattice embedding numbers. (Usually graphs represent their full weights includ- 
ing these factors.) Note that the first two graphs of the second order contribution 
also occur in a usual LCE with frozen [/-dynamics, the next two differ by an addi- 
tional dashed 2-line and the last one becomes even disconnected without the dashed 
line. 

As usual, graphical expansions for correlation functions, in particular suscepti- 
bilities, are generated from W(H, I, v) by taking derivatives w.r.t. the external fields 
H and /. Graphically this amounts to attaching external 0-lines and [/-lines with 



'1 endpoint) attached to vertices, e.g. 



(no endpoint) attached to i^-lines, e.g. 




(18) 



In passing we remark that the conventional LCE is included as a special case of the 
DLCE, if the [/-dynamics is "frozen" to some value Uq ^ 0, so that 



W\I) 
dW^I) 
dl 

dl n 



-S^UJ+IUo, 



Un 



for all n > 1, 



(19) 



i.e., no n-lines do occur with n > 1. In this case it becomes redundant to attach 
dashed lines to bare lines. As mentioned above, in an LCE only the first three 
contributions would be left in Eq. (|T 
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3 Graphical expansion 



3.1 Multiple-line graph theory 

The definition of a multiple-line graph as it will be given here is adapted to the 
computation of susceptibilities and the free energy, where the sum is taken over 
all possible locations of the fields. The definition easily generalizes to correlation 
functions. 

For details of the standard definiton of graphs in the framework of linked cluster 
expansions and related notions we refer e.g. to 0, [5]. Here, for convenience, we briefly 
recall the very definition of a graph to point out the new properties of multiple-line 
graphs as defined below in this section. 

A (standard LCE) graph or diagram is a structure 

f = (£ r ,£ r ,£r,$r), (20) 

where Cr and Br ^ are disjoint sets of internal lines and vertices of T, respectively. 
Er is a map 

Er-.Bv -> {0,1,2,...}, 

v -> E r (v) (21) 

that assigns to every vertex v the number of external lines Er(v) attached to it. 
Finally, <£>r is the incidence relation that assigns internal lines to their two endpoints. 
A multiple-line graph or multiple-line diagram is a structure 

r = {C r ,Mr,B r ,E« , \E$ n ,$ r ,*r)- (22) 
Cr, M-r and Br are three mutually disjoint sets, 

Cr = set of bare internal lines of T, (23) 
Air = set of multiple lines of T, (24) 
Br = set of vertices of T. (25) 



E T is a map 



4*>:Br - {0,1,2,...}, 

v -> E T *\v) (26) 



that assigns to every vertex v the number of bare external 0-lines E r ^\v) attached 
to v . Every such 0-line represents a field <fi. The number of external 0-lines of V is 
denoted by E^r = ^2 v€Br E r ^ (v) . Similarly, E r u ^ is a map 

E p-.M r - {0,1,2,...}, 

m -> E ( T U \m) (27) 
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that assigns to every multiple line m the number of external [/-lines E T '(m) at- 
tached to m. Every such [/-line represents a field U associated with a lattice link. 
The number of external [/-lines of T is given by E r u ^ = YlmeMr E r u \m). 

Furthermore, $r and \&r are incidence relations that assign bare internal lines 
to their endpoint vertices and to their multiple lines, respectively. We treat lines as 
unoriented. The generalization to oriented lines is easily done. More precisely, let 
(Br x Br) be the set of unordered pairs of vertices (v,w) with v,w G Br, v ^ w. 
(The bar implies unordered pairs, the prime the exclusions of (v,v), v G Br-) As in 
the standard linked cluster expansion, self-lines are excluded. Every bare internal 
line is then mapped onto its pair of endpoints via 

$r : C T -> (Br x Br) '■ (28) 

We say that v and w are the endpoint vertices of Z G Cr if $r(0 = (v,w). If there 
is such an Z G Cr, v and w are called neighbours. Similarly, \I/ r is a map 

: £r -> Mr, 

I - * r (Z) (29) 

that maps every bare internal line to a multiple line. A multiple line m G .Mr 
is composed of bare internal lines Z G Cr which belong to m in the sense that 
\E'r(0 = m. Ijvirh 71 ) lii ^e total number of bare internal lines belonging to m. With 
v — Zx r ( m ) + Er U \ m )i m 1S called a z/-line. We always require that v>\. On the 
other hand, every bare internal line belongs to one and only one multiple line. For 
simplicity we often identify a 1-line with the only one bare line that belongs to it. 

Next we introduce some further notions that will be used later. External vertices 
are vertices having external 0-lines attached, 

B r , ex t = {veBr\El?\v)^0}, (30) 

whereas internal vertices do not, Br,i n t — &r \ Br^xt- Similarly, external multiple 
lines have external [/-lines attached, 

Mr, ext = {m G Mr | E r U \m) ^ 0}, (31) 

and the complement in Mr are the internal multiple lines, Mr,int = Mr \ Mr,ext- 
For every pair of vertices v, w G Br, v ^ w, let $ (v, w) be the set of lines with 
endpoint vertices v and w, and |$ (v,w)\ the number of these lines. Thus $ (v,w) 
is the set of lines v and w have in common. With E T ^\v ) denoting the number of 
external 0-lines attached to v G Br, 

tsAv) = J2 + 4%) (32) 

weB r 

is the total number of bare lines attached to v. 
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Some topological notions and global properties of graphs will be of major interest 
in the following. A central notion is the connectivity of a multiple-line graph . Recall 
that we want to consider the DLCE expansion of the free energy and of truncated 
correlation functions as an expansion in connected graphs. As indicated in section 
2, the main generalization compared to the common notion of connectivity of a 
graph which is required here is that an additional type of connectivity is provided 
by multiple-lines. To define the connectivity of a multiple-line graph T, V first is 
mapped to a (standard) LCE graph T to which the standard notion of connectivity 
applies. There are various equivalent ways to define such a map. We choose the 
following one. 

• For every multiple-line m G Mr define a new vertex w(m). Let Br = 
{w(m)\m G Mr} and define B = B r U B T as the union of the vertices of 
T and the new set of vertices originating from the multiple-lines. 

• For every bare internal line / G Cr define two new internal lines h,l 2 and 
incidence relations 

£(Zx) = (v^w^ril))), 

where v\ and v 2 are the two endpoint vertices of I 
for all I G Cr, is denoted by C. 

• Define the external incidence relations 

E : B - {0,1,2,...}, 

E(v) = Ep\v), forveBr, 
E(v) = E r U \m), for v = w(m) G Br 

Now, T is defined by 

f = (£,£,£,$)• (35) 
Having defined the standard LCE graph T for any multiple-line graph T, we call 
T multiple-line connected or just connected if T is connected (in the usual sense). In 
Fig. 4 we have given two examples for a connected (upper graph) and a disconnected 
(lower graph) multiple-line graph . 

The next important notion is the topological equivalence of two multiple-line 
graphs . Two multiple-line graphs 

r, = (d,Mi, Bi,E^\EP,^i, i = 1,2 (36) 

are called (topologically) equivalent if there are three invertible maps 





- B 2 , 




02 : £1 - 


- A, 


(37) 


03 : Mi - 


- A^2, 
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(33) 

The set of all lines h,h, 



(34) 



Figure 4: Example of multiple-line connectivity. The upper multiple-line graph r\ 
is connected because the graph r\ is connected in the conventional sense. The lower 
multiple-line graph T 2 is disconnected because T 2 is so. 



such that 

$ 2 O 2 = (f) l O $1, 

$2° 02 = 03°^1, (38) 

and 

£f } o0 3 = ^r. (39) 

Here o means decomposition of maps, and 



(f> 1 :B 1 xB 1 ^ B 2 x B 2 

4i(v,w) = (Mv),Hw)). (40) 

A symmetry of a multiple-line graph T = (£, A4, B, E^ u \ $, \l/) is a triple 
of maps 0i : B — > £>, 2 : £ — > £ and 3 : .M ^ A"! such that 

$ o 02 = 0i o <E>, 

\I/O02 = 03 O V[>, (41) 



and 



o 0i = 

EMofa = E (u \ (42) 



The number of these maps is called the symmetry number of T. 
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We denote by Gei,e 2 {L) the set of equivalence classes of connected multiple- 
line graphs with L bare internal lines, E\ external 0-lines and E 2 external [/-lines. 
Furthermore we set 

G Ei ,e 2 ■= |J Q Ei ,e 2 (L). (43) 

L>0 

A multiple line graph T does not need to have a vertex. If Br = 0, we have 
Cr = as well. If in addition T is connected, Air consists of only one element, 
with all external [/-lines attached to it. (We anticipate that T is 1-multiple-line 
irreducible (1MLI) by definition. For the definition of 1MLI cf. section 4 below.) 
The only graph of Qq,e(L = 0) is given by 




It represents the leading term of the susceptibility 

^ = VB E <U(k)---U(l E )> c (45) 

Zi,...,Zb£Ai 

and is given by d E W 1 ( y I)/dI E \ I _ Q . The index c in ((45p stands for truncated (con- 
nected) correlation. 

By removal of a v-line m G Air we mean that m is dropped together with all bare 
internal lines and all external [/-lines that belong to m. This notion is explained in 
Fig. 5a. (It is used in section 4 for 1-lines to define 1-particle irreducible (1PI) and 
1-line irreducible (1LI) multiple- line graphs .) 

On the other hand, by decomposition of a z/-line m G Air we mean that m is 
dropped together with the external [/-lines of m, but all bare internal lines that 
belong to m are kept in the graph, being identified now with 1-lines. This notion 
will be used below to define 1MLI and renormalized multiple-line moments. It is 
illustrated in Fig. 5b. 

Similarly, decomposition of a vertex v G Br means to remove the vertex v and 
to attach the free end of every line that entered v before to a new vertex, a sepa- 
rate one for each line. This notion is used to define 1- vertex-irreducible (1VI) and 
renormalized vertex moments for multiple-line graphs . For an example see Fig. 5c. 



3.2 Susceptibilities and weights 

In the last section we have defined multiple-line graphs and the notions of connec- 
tivity and equivalence of such graphs. The definition is chosen in such a way that 
the series expansions of the free energy and of truncated correlation functions are 
obtained as a sum over equivalence classes of connected multiple-line graphs . The 
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Figure 5: Removal (a) and decomposition (b) of a 5-line. Decomposition (c) of a 
vertex. 



number L of bare internal lines of a multiple-line graph T counts the order in the 
expansion parameter v(x,y) to which T contributes. If v(x,y) is of the form 

v(x,y) = 2K J2 <U ( 46 ) 

with Af(x) any finite x-dependent set of lattice sites, the contribution of T is a 
multiple of (2K) L . Often used special cases are the nearest neighbour interactions 

D-l 

v(x,y) = 2K^2(S x ,y +Ji + S x> y- fi ) (47) 
and the uniform interaction 

v(x,y) = 2K (1-S x , y ), (48) 
which is used in models of spin glasses and partially annealed neural networks. 
Susceptibilities of the and U fields will be represented as 

Xe u e 2 = yjj Yl <<P{xi)---<p{x E ,)U{l 1 )---U{l E2 )> c 

xi,...,rr Bl eA ^,...,1^ eAi 

- 1 d El+E2 W(H, J, • 



VD x 1 t%Ao, , ,a dH(x 1 )...dH(x El )dI(l 1 )---dI(l E2 ) 



H=I=0 



= D 2ir ) L E w ( r ) ( 49 ) 

£>o reg Bl ,E 2 (L) 

with lattice volume V 7 and dimension Similar representations hold for higher 
moments /i. 
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The weight w(T) of a multiple-line graph T G Gei,e 2 {L) is given as the product 
of the following factors 



for every vertex v G Br a factor 

U " = l^ipHw (50) 
where n = tg r (f ) is the total number of bare lines attached to v. 

for every multiple line m G M.y a factor 

< - (^)w 

where 1/ = lMr( m ) + ^r^( m )j that is ^ is a z/-line, 
a factor l/5r, where Sr is the topological symmetry number of V, 
a factor counting the permutation symmetry of external 0-lines, 

#1 

T . (52) 

a factor counting the permutation symmetry of external ?7-lines, 

iu^w ' (53) 

the lattice embedding number of T, which is the number of ways T can be 
embedded on a lattice of given geometry, e. g. on a hypercubic lattice. To this 
end, the vertices of T (if any) are placed onto lattice sites. One arbitrary vertex 
is placed at a fixed lattice site, in order to account for the volume factor 1 /V 
in (|49D . A priori there is no exclusion principle. This means that any number 
of vertices can be placed at the same lattice site. (This is sometimes called free 
embedding.) Two restrictions apply to the embeddings. The first constraint 
results from the fact that a bare internal line represents a hopping propagator 
v(x, y), with lattice sites x and y at which the two endpoint vertices of the line 
are placed at. A reasonable computation of the embedding number takes into 
account the particular form of v(x, y) from the very beginning. The second 
constraint is that bare lines of the same multiple-line have to be mapped on 
the same pair of sites. 

For example, if v(x,y) is the nearest neighbour interaction (|47|), two vertices 
which have at least one line in common are to be placed at nearest neighbour 
lattice sites. On the other hand, a propagator v(x, y) of the form (]48f) implies 
a rather weak constraint in that x and y must be different, but otherwise can 
be freely placed over the lattice. 
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We remark that in case of a non-trivial internal symmetry (such as considered 
in section 7) the expressions of Eq.s fl49|)-(|5lD must be modified appropriately. In 
particular, the weight (|51"| ) of a multiple-line does no longer take such a simple form. 



4 Renormalization 

Truncated correlation functions, susceptibilities and other moments are obtained as 
sums over multiple-line graphs that are connected. Their number rapidly grows with 
increasing order, that is with increasing number of bare internal lines. The procedure 
of "renormalization" means that the connected moments are represented in terms 
of reduced ones. The reduced moments are obtained by summation over multiple- 
line graph classes which are more restricted than just by their property of being 
connected. Of course the number of graphs of such classes is smaller. Only the 
most restricted multiple-line graph classes must be constructed. The subsequent 
steps towards the moment computation are most conveniently done by operating 
analytically with the reduced moments. In particular, it is no longer necessary 
to generate all connected and the corresponding intermediate multiple-line graph 
classes. 

A connected multiple-line graph T is called 1-particle irreducible (1PI) if it sat- 
isfies the following condition. Remove an arbitrary 1-line of T. There is at most one 
connected component left that has external lines attached. (This notion is the same 
as the one used in the context of Feynman graphs.) On the other hand, if in addition 
the remaining graph is still connected, then V is called 1-line irreducible (ILI). In 
many cases it is sufficient to use only the second notion. It is for instance sufficient 
that all vertices are constrained to have only an even number of lines attached, or 
more generally, if graphs and subgraphs with one external line are forbidden. For 
notational simplicity we assume in the following that this is the case and henceforth 
refer only to the notion ILI 4 . The generalization to the case in which ILI and 1PI 
graphs must be distinguished goes along the same lines as for LCEs, which was 
discussed in [ |15|] . 

By G^e^L) we denote the subset of multiple-line graphs T e Q El E2 (L) that 
are ILL lLI-susceptibilities are defined as series in the hopping parameter similarly 
as in P9| ) by restricting the summation to ILI graphs, 

Susceptibilities are easily obtained in a closed form in terms of lLI-susceptibilities 
X 1LI ■ L can be shown that the x 1LI s can t> e obtained by an appropriate Legendre 
transform. For instance 

ILI 

A2,0 

X2 '° R«' 



J In ref. [|[ [| the term 1PI was used instead. 
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where v(k) is the Fourier transform of the hopping propagator v(x,y), 

v(x,y) = J\^n^ Hx ~ y) m- (56) 

In LCEs the second important resummation comes from so called vertex renor- 
malizations. This means partial resummation of graphs with specific properties such 
as having one external vertex only. These sums then are considered as " renormalized 
vertices" replacing the vertices of graphs with complementary properties. The pro- 
cedure naturally leads to the notion of 1-vertex irreducibility (1VI) and renormalized 
moments. 

In DLCE we follow this procedure. The very definition of 1VI has to be modified 
slightly for multiple-line graphs because of the enhanced connectivity properties due 
to multiple-lines. In addition, as a natural generalization, we supplement vertex 
renormalization by multiple-line renormalization. 

A multiple-line graph V is called 1-vertex irreducible (1VI) if it satisfies the 
following condition. Decompose an arbitrary vertex v £ Br- Every connected 
component of the remaining graph has then at least one external line attached. It 
can be a 0-line or a U -line. We write 

G^kiL) = {r g g^ E2 (L) | r is ivi} (57) 

for the set of equivalence classes of graphs that are both 1LI and IVI, with E\ 
external 0- lines, E 2 external ?7-lines and L bare internal lines. 

The renormalized vertex moment graphs are 1LI graphs that have precisely one 
external vertex and no external multiple line, 

Q k (L) = {r G Qlfl{L) | there is v G B T with E^\v) = k}. (58) 

A multiple-line graph T is called 1-multiple-line irreducible (1MLI) if it satisfies 
the following criterion. Decompose an arbitrary multiple-line m G Air- Every 
remaining connected component has then at least one external line attached. It can 
be a 0-line or a [/-line. We write 

QZeI(L) = {V G g^ E2 (L) I r is 1MLI}. (59) 

The renormalized multiple-line moment graphs are graphs that are 1LI and have 
precisely one external multiple-line, but no external vertex, 

R k (L) = {re Q] L k \L) | there is m G M r with E T U \m) = k}. (60) 
The equivalence classes of graphs that are both IVI and 1MLI are denoted by 
S Ei ,e 2 (L) = CkW^C^)- (61) 
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With the renormalized moment graphs as defined above, the lLI-susceptibilities 
are now obtained in the form 

X^ 2 = E( 2K ) L E ^ r )- (62) 

The weights w(T) are given as a product of factors as described in the last subsection, 
with the following two exceptions. 

• The vertex coupling constants v ° c are replaced by the renormalized vertex 
moments 

V°n - < = E( 2 ^) L E ( 63 ) 

L>0 veQ n {L) 

• The multiple line coupling constants m* c are replaced by the renormalized 
multiple line moments 

< - < = E( 2 ^) L E ™( r )- ( 64 ) 

L>0 T&R V {L) 



In the series representations above, the w(T) are computed according to the rules 
of subsection 3.2. 

In |TJ| we have described an algorithmic construction of graphs that is the first 
step for an automatic computer aided generation. In exceptional cases DLCEs can 
be summed up in a closed form. Otherwise, a computer aided generation of graphs is 
unavoidable, because the proliferation of DLCE graphs is pronounced even compared 
to LCE graphs. For LCE graphs we know that billions of graphs must be included 
to obtain the critical exponents with the accuracy of some per mil. 



5 Applications to spin glasses 

In this section we consider applications of DLCEs to disordered systems, in particular 
to spin glasses with '"slow"' interactions coupled to '"fast"' spins. The interactions 
J are assumed to be in equilibrium with a thermal heat bath of inverse temperature 
(3\ while the spins a are equilibrated according to a second inverse temperature 
(3. Both systems need not be mutually in equilibrium. Let Zp(J) be the partition 
function that describes the equilibrium distribution of the spins for given Js, 

,J(ij)(Ti<7j). (65) 

{cn=±l} i<j 

The sum runs over pairs that need not be restricted to nearest neighbours only. 
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We further assume that the dynamics of the time evolution of the slow interac- 
tions J is governed by a Langevin equation 

with 

1 1 - 

H(J) = --\nZ p (J) + -iiN J hv ( 67 ) 

P i<j=l 

Here Zp(J) is given by (|65|) , iV is the total number of spins, \i is a positive constant 
and 77jj is a stochastic gaussian white noise of zero mean and correlation 

< Vv(t)Vki(t') > = jAiMkiAt - t'). (68) 

Such a Langevin equation for the Js can be derived from an ansatz which is mo- 
tivated by neural networks |Il|- ||13|| . Moreover, since the time evolution of the Js 



is determined by a dissipative Langevin equation, the equilibrium distribution of 
the slow variables is again a Boltzmann distribution, governed now by the second 
temperature 

poo n 

Z' v = J\f J J] dJ (i>j) exp(-/3'ft(J)), (69) 

i<j=i 

with M some normalization that will be specified below. The effective Hamiltonian 
TC of J is given by (p7|). 

It is these equilibrium aspects of coupled systems of fast spins and slow interac- 
tions that we can treat analytically with DLCEs, as we will show below. 

Let us first rewrite ZL in the form 

Z '<# = r fi \\fW dJ{i > j n -^^ N zZ J h)) z m x 

00 i<j=l V / i<j 

= [[MJ)% (70) 

where we have introduced Q = (3 fx and real x = (3' / j3 as the ratio of two tem- 
peratures. The normalization has been chosen such that [[1]] = 1. In the limit of 
x — > for fixed (3, i.e. (3' — ► 0, we have a quenched system. The J only feel the 
infinitely high temperature but are decoupled from the spins. The time scale 

of fluctuations of the spins is assumed to be so short that the J are only sensitive to 
averages of the a. Vice versa, the spin dynamics does depend on the Js. Therefore 
the quantity of physical interest is not 

In 



■°°i<j=l V / Wi} i<j=l i<j 

(71) 
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where fluctuations of the as do influence the Js, but 

N / H^T7 \ N 



J II ( V27 rfj ^)) ln [E exp (^ E J {ij) a i a i) - ex v(-lQ N ^2 J h 

(72) 



x i<j=l \ ' J U {cn} i<j=l ' i<j 



or, in a shorthand notation, 



J VJ ln^(J) = [[fnZ^(J)]]. (73) 



Usually one rewrites 

z p {jy - 1 



/ 



VJ lnZfl(J) = I VJ lim 
J 



x 



fa/ PJ MSL-L = lim Mi + (gW11-i)} 



InZ' 

= lim ^. (74) 

For the second equality sign one has assumed that J VJ commutes with lim^o, in 
the third one that lim :r ^ [[Z /3 ( J) x ]} = 1. So far, x as the ratio of two temperatures 
is real. Rewriting the left hand side of ( [74]) according to the right hand side is 
called the replica trick 17]. The uncontrolled approximation that usually enters 
the replica trick is that now the right hand side is evaluated for positive integer 
x = n and extrapolated to n — 0. Clearly a function that is known only for positive 
integer n does not have a unique extrapolation to n = without further assumptions. 
Nevertheless, this approximation is made, because it is rather convenient. For integer 
n, Zp(J) n is the partition function of an n times replicated system of which the 
logarithm is taken after the integration over the Js. It is seen as follows. We rewrite 

n N 

MJ) n = E ex p(^E E J w)^ a M 0) )> ( 75 ) 

with a = 1, . . . , n labelling the replicated spin variables, so that 



/oo ly 
J] dJ [itj y E exp(-5(J,a^)), 
°°i<i=i K< a >=±i} 

n iV 1 

5(j^ (o) ) = -/?E E %^ fl M fl) + 2^E4)' w 

a=l i<j'=l i<j 

Linear terms in cr and J may be included according to 

n 2V N 

S lin = -^EE^ a) + c E J (^) (77) 

a=l i=l i<j'=l 
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with constant external fields h and c. 

Apparently, because of integer n, Z' n/3 has the form of models to which DLCE 
applies, with a hopping term 

n N 

S hop (J,° (a) ) ="^EE %^N a) , (78) 

a=l i<j=l 

a single link action 

S\J {iJ) ) = cJ {iij) + l -QN J^, (79) 

and a single site action 

n 

5°(^ (a) ) = -hJ2^- (8°) 

a=l 

Depending on n we distinguish the following cases. 



n — 1. First we note that for n = 1 we can directly apply DLCE to lnZ^, =l3 
and to derived quantities to obtain their series expansions in (3. But from a 
physical point of view, in a disordered system one is not interested in n = 1, 
because n = 1 corresponds to the completely annealed situation, in which the 
fast spins and the slow interactions are in mutual equilibrium. (In contrast, 
in particle physics one is interested in the n = 1 case, cf. our applications 
of DLCEs in the framework of variational cumulant expansions of the SU(2) 
Higgs model fT4|| .) 

n > 1, integer. Again we apply DLCE to In Zg /=nf3 , but have to account for the 
permutation symmetry between the replicas. Formally, the replica symmetry 
plays a role similar to an internal symmetry, e.g. an O(N) symmetry in a 
scalar Higgs model. DLCEs with nontrivial internal symmetries have been 



discussed in connection with the SU(2) Higgs model Thus we can study 
'"unquenched"' equilibrium aspects of systems with two temperatures and 
compare the results from DLCEs adapted to "internal" replica symmetry with 



Monte Carlo simulations for the same n fllq, 19 



n = 0, the quenched limit. As we will show in the next section, in order 
to discuss the x — > limit, we need not refer to n times replicated systems 
Z p (J) n characterized by @-(|0D, but just to Z p (J) given by ©-® with 



n = 1. By means of special DLCEs , so-called quenched DLCEs (QDLCEs), 
we directly calculate the left hand side of fl74l). Therefore, setting n = 1 in 
([7|)-(^(]) in QDLCEs does not imply the completely annealed case, because 
we first take the logarithm of Zp(J) and then average over the Js. 
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5.1 Avoiding the replica trick 

First we adapt the notation to section 2 to include more general cases. Ao denotes 
the support of the spins, that is the set of lattice sites, with V = \Aq\ denoting their 
total number. Ai C Ai are the pairs of sites whose spins interact. In accordance 
with (ffOD , we write for the normalized link- average of a function f(J) 

mm = j vjf(j) (si) 



with 



vj = Hdfi(j(i)) , 



leAi 



/oo 
dfi(J) = 1. (82) 
-oo 

It is convenient to introduce the single link expectation values 

< g(J) >i = J dfi(J) g{J) (83) 

and the generating function W X {I) by 

exp W\I) = < exp (I J) >i . (84) 

The way in which the replica trick can be avoided is examplified for the free 
energy density W sp /V of the spin system averaged over the link couplings. The 
partition function of the spin system for a given distribution of the link interactions 
J(x, y) is given by 

expWVJ) = K P J Va exp (—S sp (a, J)), (85) 
where W sp (0) = and 

S sp (a,J) = - - ^ v(x,y)a(x)a(y)J(x,y), 

x,y£A 

Vo = Y[ da(x) -exp(-S (a(x))). (86) 

xeA 

Without loss of generality we identify the support of the interaction v(x, y) = v(y, x) 
with the set A 1 of lattice sites where VJ is supported, 



Ai = {I = (x, y) G A x A | v{x, y) ^ 0}. (87) 

For simplicity we assume v(x,y) to be of the form fl46|), so that K is a measure of 
the strength of the interactions v(x, y). 
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The free energy density of the spin system allows for a series expansion in the 
standard LCE sense, with the link field J (I) playing the role of a "background field", 



±w sp (j) = E(2^) L E ^OV). (88) 

l>o reg s p (L) 

Here Q S £{V) (with E = 0) denotes the set of equivalence classes of connected LCE 
graphs with E external lines and L internal lines. The spin- weights w sp (T, J) are of 
the form 

w sp (T,J) = R sp (T) U J ( l ) m(l) - ( 89 ) 

The sum is taken over all non- vanishing lattice embeddings of the graph T. It runs 
over all maps of internal lines of the graph T to pairs of lattice sites of Ai that are 
consistent with the graph topology in the sense discussed in section 3. For every 
/ G Ai, m(l) denotes the number of lines of T that are mapped onto the link / by 
the embedding. All other factors that contribute to the weight are collected in the 
prefactor R sp (T), including the inverse topological symmetry number of T. 

Next we want to express [[W sp ( J)}} as a series in K by means of DLCE. Toward 
this end we set f(J) = W sp (J) and insert the series (]8§D with (|39| ) into (|8lD . At this 
stage we are not concerned with question of (uniform or dominated) convergence 
and obtain 



[[^W sp (J)\] = E( 2 ^) L E hjw sp {T,J) 



= E( 2ir ) L E RSP ( r ) E n < J (o m(0 >i • ( 9 °) 

The next step is to express the full single link expectation values in terms of the 
connected ones. They are related by 

< J m >i = e n < j|p >s. ( 9i ) 

nep(m) Pen 

where V(m) denotes the set of all partitions of m = {1, . . . , m} into non-empty, 
mutually disjoint subsets of m. \P\ is the number of elements of the set P. The 
relation fl9l| ) is equivalent to the partition of all lines of T that are mapped to the 
same lattice link into multiple lines, with every multiple line contributing a factor 



dI\P\ 



lc 

PI 

/=0 



™\p\- (92) 
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Using we rewrite @0|) as 

[[^(J)]] = £(21C) £ £ i^(r) £ (iRpi) (EII 1 

l>o reg^(L) neP(£ r ) \Pen / yn-^Xi /eAi 

(93) 

The last summation in ( |93| ) is over all maps Cr — > Ai of the lines of T to the 
lattice links of Ai subject to the constraint that all lines that belong to the same 
multiple-line corresponding to some P G II are mapped onto the same lattice link. 

Finally we rewrite (|93|) as a sum over multiple-line graphs . To this end, we 
first observe that for every T G Qq P (L), every partition II G V{Cr) of the lines 
of r into multiple-lines generates a multiple-line graph A = (r, II) in the obvious 
way. Let us denote by Go,o(L) the subset of multiple-line graphs of Go,o(L) that 
stay connected after decomposition of all multiple lines. (These are the multiple- 
-line graphs which stay connected in the usual graph theoretical sense, when the 
dashed lines are omitted.) For every A G Go,o(L) there is a unique T(A) G Gq P {L) 
and at least one II G V(Cr(A)) such that (r(A),II) = A. Let be the (uniquely 
determined) number of partitions II G V(Cr{A)) with (r(A),IT) = A, and 11(A) 
such an arbitrary partition. Eq. (|93|) then becomes 



[[^W sp (J)]] = J2m) L E n A R^(T(A)) ( J] m $\ ( ^_ 
l>o Aeg 00 (L) \Pen(A) / \n(A)— Aj 

(94) 

The last bracket of ( p4] ) is the lattice embedding factor of the multiple-line graph 
A. The second bracket from the right does not depend on the choice of 11(A) and is 
the product of the multiple-line coupling constants as defined in section 3. Finally, 
?7.A-R' sp (r(A)) is precisely the remaining part of the weight of A that was described in 
detail in section 3, endowed with the correct inverse topological symmetry number 
of the multiple-line graph A (because of the factor n&). 

In summary, we obtain the series expansion of the link-averaged free energy 
density in terms of DLCE graphs, 

4 W -^ J ^ = E( 2 ^) L E W (A). (95) 
L>o Aeg 0fi (L) 

The weight it? (A) of a multiple- line graph A is defined and computed according to 
the rules given in section 3. 

Eq. ( p5|) is the series representation of the link-averaged free energy density of 
the spin system, i.e. the free energy density of the n = replica system, in terms of 
DLCE graphs. It looks much like the series representation of the 1-replica system, 
which is given by 

^W^ repl = i ]n[[exvW sp (J)}} = E( 2 ^) L E w ^ ( 96 ) 

L>o AeS ,o(£) 
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according to the discussion of section 2. We recall that Go,o(L) is the set of DLCE 
vacuum graphs with L bare lines that are connected in the generalized DLCE sense. 
Comparing ( PS|) and (|96|), the transition from n = 1 to n = replicas is achieved by 
keeping only the subset Q ^(L) C £ (L) of multiple-line graphs that are connected 
in the original (LCE) sense. 

We emphasize that the restriction of DLCEs to QDLCEs is not an ad hoc (or 
intuitively motivated) assumption but a derived consequence of the fact that the 
logarithm is taken before the integration J T>J . This procedure accounts for all 
graphs that contribute to a given order in K. Thus we do have to truncate the 
series unless the series can be completely summed up, as it happens in exceptional 
cases. 

We expect that the series (^) are convergent for a large class of interactions 
5' 1 (J) and v(x,y) if the coupling constant K is sufficiently small. For special inter- 
actions most of the multiple-line graphs yield vanishing contributions so that we can 
further restrict the sum to a subset of Qq^{L). An example is given by the mean 
field type of interaction of the Sherrington-Kirkpatrick model, cf. Sect. 5.2. 

5.2 Applications of QDLCEs 

In the following we list some examples for systems of which we can study the phase 
structure by means of QDLCEs. Their actions are special cases of ([T^)-([5D|) with 
n = 1 (as explained above) and the following choice of variables. 

• Infinite range models. Choose J(jj) G IR as before, a-i G ±1, i, j G {1, . . . , iV}, 

K(l-6 x>y ), 

Ai is the set of all pairs of sites, (97) 

N - J 2 . 
2 

Now the sum over the sites in ( |78| ) runs over arbitrary pairs (i, j), i < j, and 
we obtain the infinite range Sherrington-Kirkpatrick model. For infinte range 
and in the thermodynamic limit (N — > oo), the phase structure can be solved 
by replica mean field theory, cf. e.g. [pl|. For QDLCEs the infinite range and 
N —>■ oo limits imply that only tree graphs of 2-lines contribute to the series 
of the free energy density, such as 



(98) 



v(x,y) = 
Ai = 
S\J) = 
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The reason for that is that each 2-line gets a factor of 1/N from the Si- 
part of the action, but each vertex gets a factor (N — 1) from the embedding 
onto a lattice Ao- The contribution of every tree graph to the free energy is 
proportional to N. If a chain of 2-lines connecting the vertices gets closed, 
forming a loop, there is one (N — 1) less in the total embedding factor. Thus 
the contribution is suppressed by 1/ (N — 1) for every loop and vanishes in the 
thermodynamic limit. Because of the simple tree structure there is a chance 
for summing up the series. This is currently under investigation. 

• Finite range connectivity. The sum Yli<j °f the spins is now restricted to 
next-neighbours or, more generally, to a finite number of pairs. Rather than 
specifying Si(J(jj)) of (|79|), it is sufficient for DLCEs and QDLCEs to choose 
exp(-Si(J (i)j) )). Let 

exp{-Si{J)) = (l-p)5{J)+ P 5(J-l) (99) 

with p G [0, 1]. The variables Ju^ G {0, 1} can then be interpreted as occu- 
pation numbers of the bonds. Furthermore, if we choose <7j G {±1} we obtain 

a 

— bond-diluted Ising model. 
Choosing G Z g we obtain a 

— bond-diluted g-state Potts model. 
If <7j G S q , we obtain a 

— bond-diluted Heisenberg model. 

QDLCEs provide a systematic analytic expansion for disordered systems with 
bond dilution in a quenched limit. Coming from the high temperature (small /3) 
region one can study the phase structure as a function of the degree of dilution. 
Work in this direction is in progress. 

6 Summary and Conclusions 

In this paper we have introduced a new expansion scheme for 3-point interactions 
or, more precisely, for point-link-point interactions. This scheme generalizes linked 
cluster expansions for 2-point interactions by including hopping parameter terms 
endowed with their own dynamics. In chapters 3-4 we have developed a multiple- 
line graph theory with an additional new type of multiple-line connectivity. We 
have introduced appropriate equivalence classes of graphs and discussed the issue of 
renormalization. These notions are required for an algorithmic generation of graphs. 
Because of the fast proliferation of graphs already at low orders in the expansion, a 
computer aided implementation becomes unavoidable, if one is interested in higher 
orders of the expansion than we have computed so far. 
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In Sect. 5 we have shown how to avoid the replica trick for calculating the free 
energy of disordered systems in the quenched limit. DLCEs are a systematic expan- 
sion method to study the phase structure of disordered systems. It is systematic in 
the sense that we do not restrict the expansion to certain subclasses of graphs that 
can be summed up, but we identify and keep all graphs that contribute to a given 
order in the expansion parameter. DLCEs provide an analytic tool for studying 
systems in situations in which it has been impossible so far. 
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